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ABSTRACT 

Molecular hydrogen H2 is a key requirement for star-formation, frequently observed along sightlines in the 
Galaxy and to the Magellanic Clouds, but notoriously hard to detect beyond z ~ 0. Because H2 has no dipole 
moment, and its absorption signatures, the Lyman-Werner bands, are located bluewards of the Lymana tran- 
sition, H2 is most directly probed at z > 2 only through high-resolution spectroscopy of damped Lymana 
(DLA) systems. In the DLAs associated with long y-ray bursts (GRBs), which are tightly linked to vig- 
orous star-formation, H2 has, however, remained largely elusive. Molecular gas in GRB-DLAs has been 
unequivocally detected only towards GRB 080607 and possibly towards GRB 060206. Here, we present 
the discovery of H2-rich gas, including the presence of vibrationally-excited H^ in the optical spectrum of 
the afterglow of GRB 1208 15A at z - 2.36 obtained with X-shooter at the VLT. Simultaneous photometric 
broad-band data from GROND and X-ray observations by Swift/XRT place further constraints on the amount 
and nature of dust along the sightline. The galactic environment of GRB 1208 15A is characterized by a 
strong DLA with log(A/(Hi)/crrr 2 ) = 21.95 ± 0.10, prominent H2 absorption in the Lyman-Werner bands 
(log(/V(H 2 )/cnr 2 ) = 20.53 ± 0.04) and thus a molecular gas fraction log/(H 2 ) = -1.14 ± 0.10. The dis- 
tance d between the absorbing neutral gas and GRB 1208 15 A is constrained via photo-excitation modeling 
of fine-structure and meta-stable transitions of Fe II and Ni II to d — 0.5 ±0.1 kpc. The DLA metallicity 
([Zn/H] = -1.15 ± 0.11), visual extinction (A v < 0.15 mag) and dust depletion ([Zn/Fe] = 1.01 ± 0.09) are 
intermediate between the values of well-studied, ^-deficient GRB-DLAs observed at high spectral resolution, 
and the approximately solar metallicity, highly-obscured and H2-rich GRB 080607 sightline. With respect to 
/V(Hi), metallicity, as well as dust-extinction and depletion, GRB 1208 15 A is fairly representative of the av- 
erage properties of GRB-DLAs. This demonstrates that molecular hydrogen is present in at least a fraction of 
the more typical GRB-DLAs, and H2 and H* 2 are probably more wide-spread among GRB-selected systems 
than the few examples of previous detections would suggest. Because H* 2 transitions are located redwards of 
the Lymana absorption, H* 2 opens a second route for positive searches for molecular signatures also in GRB 
afterglows at lower redshifts and observed at lower spectral resolution. Further detections of molecular gas 
in GRB-DLAs would allow statistical studies, and, coupled with host follow-up and sub-mm spectroscopy, 
provide unprecedented insights into the process and conditions of star-formation at high redshift. 
Subject headings: Gamma-ray burst: individual: GRB 1208 15A — galaxies: high-redshift — ISM: molecules 
— dust, extinction 
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1. INTRODUCTION 



The formation of stars is tightly co rrelated with the pres- 



ence of molecular h ydrogen H2 (e. g., Shu e t al.|1987 Blitz 
1993]|McKee & Ostriker|2007| |Bigiel et al. 2008|TmtheTo^ 
cal Universe, for example, star-formation proceeds quite ex- 
clusively in cold molecular clouds, and m olecular gas is ubiq- 
uitous in the inte rior of the Galaxy (e.g., |Savage et al.|1977[ 



Dame et al . 2001 ) or the Magellanic Clouds (e.g., Tumlin 
son et al. 2002). Similarly, at higher redshift, the increased 
star-formation rates of non-merging, z ~ 1 - 2 galaxies are 
likely related to an elevated ga s-mass fraction (e.g., Tacconi 
[eTaT1[20T0l |Daddi et al-pHO] ). Direct detections of H 2 be- 
yond z ~ are, however, challenging from the ground, and 
its presence and state at high redshift are most commonly in- 
ferred indirectly via tracers such as CO or HCN. The fact that 
these molecules have dipole moments makes them much more 
amenable to detection via mm/sub-mm spectroscopy, but in- 
ferences about H2 then rely on poorly-known conversion fac- 
tors. 

Directly studying the absorption signatures of molecular 
hydrogen, the Lyman-Werner bands, at high redshift requires 
a bright back-light such as quasi-stellar objects (QSOs) or 
y-ray bursts (GRBs) to enable medium-to-high resolution^] 
spectroscopy of damped Lya (DLA) sy stems bluewards of 
the Lya transition (e.g., |Ledoux et al. 2003). But even when 
probed in very high detail with QSO-DLAs, detections of 
molecular gas are scarce: H2, for example, was detected only 



along 13 out of 77 QSO sightlines in the sample of Noter 
|daeme et al.| (|2008|l, while observations of other mo lecules 
like HP (|Varshalovich et al.|2001||Malec et al |20T0| > or CO 
( [Srianand et al.|2 008; Noterd aeme et al. 2009 2011) remain 
limited to very few systems at high redshift. 

In the DLAs associated with long GRBa^l mol ecular gas 
remained undetected for a long time as well(e.g., Tumlinson 



et al.|20O7 Ledoux et al.|2009"). Given the very high column 



densities of gas regu larly det ected in GRB afterglow spec- 
tra (e.g., |Jensen et al.pOO l ; Hjort h et al ||2003b| IVreeswijk 



et al. 2004; Jakobsson et al. 2006 ) and their direct link to star- 



formation (e.g.,jGalama et al. 1998; Hjorth et al. 2003a), this 
' isence of H2 was initially quite surprising^ 
The incidence of H2 along QSO sightlines 



Levshakov 
& Varshalovich"l[T9"8"5| [Srianand & Petitjean||T998| |Petitjean 



et al.|2000[) is positivefy correlated with metallicity and dust 
depletion (Ledoux et aT| |~2003 ; Petitjean et al. 2006; Noter- 

hotons in the rest-frame 
and the small sample 



daeme et al. 2008| l. Dust absorbs pi 
ultra-violet (UV) very efficiently^] 



2 ' High spectral resolution is required to disentangle tb features from the 
Lya forest. 

22 GRBs are typically divided into long and short events with different 
progenitor channels. In this work, we focus exclusively on long GRBs, typi- 
cally defined as having durations Tgo > 2 s (see, e.g., Kouveliotou et al. 1993 
for details). Tgo is the time during which 90% of the GRB's y-ray photons 
are emitted. 

23 Because the material probed by GRB-DLAs is in many cases located 
at distances of several hundred parsecs or more from the y-ray burst (e.g., 
IVreeswi jk et al. 2007 20lT| |D'Elia"er al. 2009 Sheffer et al. 2009 1, the after- 
glow's UV flux itself will not photo-dissociate the Hi of the GRB-DLA (e.g., 
Tumlinson et al. 2007 , Ledoux et al. 2009 1. A genuine deficit of molecu- 
lar hydrogen in GRB-DLAs as compared to QSO-DLAs for example, would 
instead indicate that H2 is absent already prior to the burst. This could be 
caused by substantial background UV-radiation fields in their host galax- 
ies, which are capa ble of dissociating hydrogen molecules efficiently (e.g., 
Whalen et al. 2008 1. 

24 In UV-steep reddening laws similar to the one observed towards the 
Small Magellanic Cloud (SMC) and often along GRB sightlines, for exam- 
ple, the ratio between the extinction in the Lyman-Werner bands and visual 



(~10— 15 events) of GRB afterglows observed with echelle 
spectrographs with a very high resolving power (R > 40 000) 
is thus restricted to the most luminous and least obscured end 
of the afterglow brightness distribution. Consequently, the 
best studied GRB-selected absorbers are typically those with 
comparatively lo w dust and metal content ( Prochaska et al. 
|2007b[ |LVElia et al.|2009l [20101 [DeCia et al.|20l2| i. Given 
their low metal abundance and dust depletion, the previous 
non-detections of 1L in GRB-DLAs observed at high spectral 
resolution are thus consistent with QSO-DLA number statis- 
tics ( [Ledoux et a l. 2009). 
The full sample of GRB afterglows, including those ob 



served with lower spectral resolution (e.g., Fynbo et al. 2009 1, 
however, is diverse and includes direct evidence for meta 



rich gas (e.g., Watson et al. 2006; Prochaska et al. 2007a. 
Eliasdottir et al. 2009; Savaglio et al. 2012). In addition sig- 
nificant dust-depletion factors (e.g., Savaglio & Fall 2004. 
Thone et al. 2013 ; Hartoog et al. 2013), large amounts of red- 



dening (e.g., |Gre iner et al. 201 1 ; Zafar et al. 2012 ) and strong 
metal lines (e.g., Savaglio et al. 2003 ; Christensen et al. 201 1 



de Ugarte Postigo et al. 2012) are regularly detected towards 
GTJBT 

A significant fraction of GRB-selected absorbers observed 
to date are thus promising sites to search for the presence of 
molecular gap] Its detection, however, has proven to be dif- 
ficult, because the two conditions for probing H2 in GRB- 
DLAs seem mutually exclusive: instrumentally, high spectral 
resolution and signal-to-noise ratio (S/N) in the rest-frame UV 
(and thus bright targets) are required to disentangle H2 transi- 
tions from the Lya forest and physically, afterglows probing 
environments with significant metallicity and dust-depletion 
(favoring more obscured, UV-faint targets). 

In fact, H2 signatures were secure ly detected only for one 
(GRB 080607, Prochaska et al. 2009), and possibly for a sec- 
ond (GRB 060206, Fynb o et al [2006||Thone et"aL[2 008 ) GRB 
sightline. GRB 080607 was an exceptional event. It had one 
of the intrinsical ly most luminous afterglows ever discovered 
(Perle y et al.|2011 1, and was observed spectroscopically only 
minutes after the initial trigger when the transient was still 
bright ( |Prochaska ef al. 2008b). The GRB 080607 afterglow 
spec trum is cha racterized by molecular absorption of H2 and 
CO ( Prochaska et al. 2009 ) including vibrationally-excited H* 2 



r-; 1 ( Shelter et al. 



|2009p at z - 3.04. It further exhibits large 
amounts of neutral gas (log(Af(Hi)/cnr 2 ) ~ 22.7) around so- 
lar metallicity, a very high dust content of Av ~ 3.3 mag and 



evidence for a 2175 A feature ( Prochaska et al.|2009 Perley 
|et al |2011 ) . 

In this work, we present the second definite detection 
of molecular hydrogen in a GRB-DLA. GRB 120815A at 
z - 2.36 was observed with X-shooter in the spectral range 
from the UV-cutoff to 2.5 yum at medium resolving power 
(R ~ 5000 - 10000), providing both the S/N and spectral 
sampling necessary to investigate the physical conditions in 
the GRB-DLA. Our spectroscopic data are complemented by 
simultaneous optical/near-infrared (NIR) photometry and X- 
ray observations, providing strong constraints on the redden- 
ing along the sightline. 

Throughout the paper we quote errors at the Icr confi- 
dence level, report magnitudes in the AB system, and provide 

extinction A lmo ^/Ay exceeds a factor of 6. Even small values of Av thus 
provide considerable suppression of the observed flux in the rest-frame UV. 

25 Most GRB-DLAs indeed show log(W(Fe) dus 7cnr 2 ) > 14.7 (De Cia et 
al., in preparation), which is the critical value above which molecules have 
been observed along QSO sightlines (Noterdaeme et al. 2008 1. 



wavelengths and redshifts in a vacuum, heliocentric refe rence 
frame. Th e atomic data used in this work are those of lMor-l 



ton 



(2003 ) with few exceptions^] Solar reference abundances 
were taken from Table 1 of Asplund et al. (2009) and follow 
the suggestions of Lodders et al. (2009) whether abundances 
from the solar photosphere, meteorites or the average between 
the two values are used. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Swift Observations 
The Burst Alert Telescope (BAT; |Barthelmy et al.|[2005]) 



onboard the Swift satellite (Gehrels et al. 20041) triggerei 



on the long and soft GRB 120815A on^UT2-08-15 at 



4) t 
12-C 



T = 02:13:58 UTC ( [Pagani et . al.||20T2l JMarkwardt et al. 
2012} . Swift's two narrow-held X- ray an d UV/opt ical in stru- 
ments, the X-ray Telescope (XRT; Burrows et al. 2005) and 
Ultra- Violet Optical Telescope (UVOT; Roming et al. 2005 ) 
started observing the GRB field approximately 2.7 ks after 
the trigger, and dete cted the fa ding afterglow associated with 
GRB 1208 15A ( |Kennea et aLp 012 ; Holland & Paganii 20T2"l l. 
X-ray data were retrieve d "From the XRT online repository 
( [Evans et al.p007l [2009) , and analyzed in XSPEC vl2 . 7 . 8 
(Arnau3 [l996| ). Due to the faintness of the afterglow at the 
time of the observation, imaging data from UVOT have lim- 
ited S/N, are not constraining and are thus not used in the 
further analysis. 

2.2. GROND Optical/NIR Photometry 
The Gamma-Ray burst Optical/NIR Detector (GROND; 



Greiner et al. 2007, 2008) initiated automatic follow-up ob- 



servations of GRB 120815A on 2012-08-1 5 02:14:24 UTC, 
with the first data taken ~150 s after Tq (Sudilovsky et al. 



2012 1. Simultaneous imaging in four optical (g'r'i z') and 
three NIR (JHK S ) filters was performed continuously until 
05:47 UTC. Our analysis of the GROND data focuses on the 
broad-band spectral energy distribution (SED) and abs olute 
flux scale contemporaneously with the X-r ay (Section 2.1 1 
and optical/NIR spectroscopy (Section 2.3 1, as well as the 
afterglow's UV-photon y ield for the modeling of the excited 
absorption lines (Section 4.1). A detailed discussion and the- 
oretical interpretation of the afterglow multi-color light curve 
is beyond the scope of this work. 

GROND data were reduced and analys ed within 
pyraf/IRAF ( Tody 1 1993 ) in a standard manner (Kriihler et al. 
|2008| l. The transient is detected with high S/N in each filter. 
Using USNO field stars as astrometric reference, we derive a 
position for GRB 120815A of R.A. (J2000) = 18:15:49.83, 
Dec. (J2000) = -52:07:52.5 with an absolute accuracy of 
073 in each coordinate. The photometric solution w as tied 
to the magnitudes of stars from the SDSS catalog (Aihara 
et al. 2011) in g'r'i' z' observed dir ectly before and after the 
GRB field, and 2MASS field stars ( ]Skrutskie et"aLl[2006[ ) in 
JHKfi. Based on the scatter of individual calibration stars 
we estimate our absolute photometric accuracy to be 4% 
in g'r'i'z', 5% in J and H and 8% in K s , dominating the 
total p hoto met ric error. The GROND data are reported in 
Tables Al and A2 and the multi-color light curves that were 



Atomic 



Tripp 



2006 



data 
and 



for 

Ni 



Ni II(/!1370) are from | Jenkins &] 
II F9/2 as well as Fe II F9/2 from 
nttp: //www.cfa . harvard . edu/amp/ampdata/kurucz2 3/sekur . html. 
We note that we apply a similar modification to Ni II 4 Fc>/2 as in Vreeswijk 
|et al.| (2007 . In the photo-excitation modelling, we also use atomic data 



provided by 



Verner et al. 1 1999 1 and references therein. 



used a s inp ut to the photo-excitation modeling are shown in 
Figure [AT] 

2.3. X-shooter Optical/NIR Spectroscopy 

Spectroscopic observations of the GRB 1208 15 A afterglow 
in the wavelength range between 3000 and 24 800 A com- 
menced on 2012-08-15 at 03:55 UT (6.06 ks after the BAT 
trigger) with the cross-d ispersed echelle spectrograph X- 
shooter ( |Vernet et al. |20 1 1 [ > mounted at ESO's Very Large 
Telescope (VLT) UT2. They consisted of four nodded ex- 
posures in the sequence ABBA with exposure times of 600 s 
each, taken simultaneously in X-shooter's ultraviolet/blue 
(UVB), visible (VIS) and near-infrared (NIR) arms. The av- 
erage airmass was 1.37 and the mid-time of the stacked ex- 
posure is 04:17:50 UTC (7.43 ks after the BAT trigger or 37 
minutes in the rest-frame at z - 2.36). Sky conditions were 
clear and dark with a median seeing of 076 during the ob- 
servations. X-shooter spectroscopy was performed with slit- 
widths of 170, a/9 and a/9 in the UVB, VIS and NIR arm, re- 
spectively. The resolving power R - A/AA was derived from 
unsaturated, single telluric lines in the red arms, and scaled 
to the UVB arm using a standard dependence of seeing with 
wavel eng th (see Fynbo e t al.|2011| for details) and tabulated 
valuesFTof R for different X-shooter slit widths. This results 
in values of R ~ 6000, 10400, 6200 for the UVB/VIS and NIR 
arm, respectively. Because of the good seeing, this is some- 
what better than the nominal values for the given slit widths 
derived from arclamp frames (R = 5380, 8800, 5300). 

X-shooter dat a were reduced wit h the ESO/X-shooter 
pipeline vl . 5 . 8 (Goldoni et al.|2006| ), rectifying the data on 
an output grid with a dispersion of 0.15 A per pixel in the 
UVB, 0.13 A/pixel in the VIS and 0.5 A/pixel in the NIR 
arm, respectively. The wavelength-to-pixel conversion was 
chosen as a compromise between slightly oversampling X- 
shooter's spectral resolution, while at the same time limiting 
the correlation of adjacent pixels in the rectification step. The 
wavelength solution was obtained against arc-lamp frames 
in each arm, leaving residuals with a root mean square of 
0.03/0.04/0.3 A, respectively, corresponding to 3/2/5 km s~ : 
at 4000/7500/16500 A. Given the projected proximity of the 
GRB to the USNO star 0378-0977505 and to maximize S/N, 
each of the four frames in the individual arms was reduced 
separately, with the sky estimated in regions free of source 
counts. The sky region is limited because of the short slit 
length of X-shooter (11") and artifacts at the order edges, af- 
fecting the accuracy and quality of the sky-subtraction in the 
NIR arm. 

The spectra were extracted from the combined frame 
in a variance-weighted method, taking into account the 
wavelength-dependence of the centroid and width of the trace. 
To estimate the error in the combined spectrum, we used 
the noise model of the X-shooter pipeline for each of the 
four individual frames propagated through the stacking and 
extraction process. The median S/N per resolution element 
AA. is ~ 15 (4200 - 5700 A), ~ 16 (5700 - 10000 A), and 
-9/14/16/15 in Y/J/H/K. 

Flux-calibration was performed against the spectro- 
photometric standard LTT7987 observed starting on 2012-08- 
15 at 05:00 UTC. The atmospheric dispersion correctors of 
X-shooter in the UVB and VIS arm were switched off af- 
ter a technical problem. Because of the substantial airmass 

27 http: //www. eso.org/sci/facilities/par anal/instruments/ 
xshooter/inst/ 
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Figure 1. NIR-to-X-ray spectral energy distribution and model for the after- 
glow of GRB 120815A at 7.3 ks after the trigger. Solid lines show the model 
including gas and dust absorption, while the dashed line illustrates the under- 
lying synchrotron continuum emission. X-ray data have been binned to yield 
a S/N of at least 8 to enhance clarity. The g'-band photometry is not fitted, 
because this filter extends bluewards of the Lya transition. 

during our observations we expect some chromatic slit-loss 
in the data t hat w e aim to correct via our afterglow model 
(see Section |3.1). The flux-calibrated X-shooter spectrum 
low-order polynomial (iteratively exclud- 



3.1 
was fitted with , 



ing the absorption lines), which yields the continuum level 
of the data. The ratio between continuum fit and afterglow 
model then provides wavelength-dependent matching factors 
over the whole spectral range of X-shooter, and thus an accu- 
rately flux-calibrated spectrum. 

3. RESULTS 

3.1. Broad-band Spectral Energy Distribution 

We fitted the broad-band afterglow data from XRT and 
GROND under the assumption that the underlying afterglow 
flux is well represented by synchrotron emission and the af- 
terglow is in the slow-cooling regime (see 



e.g., 



Granot & 



Sari||2002] > with the optical/X-ray energy range being above 
the afterglow's characteristic synchrotron frequency. Single 
and broken power-law models were used, where the soft X- 
ra y c olumn density is attributed to two absorbers with so- 
laiF] metallicity : one Galactic he ld fixed at JV(H) X = 8.6 x 



lCFcnr 2 (Kal berla et al.||2005) , one intrinsic to the GRB 
host. The reddening in the optical/NIR wavelength range is 
modeled by extinction laws from the Magellanic Clouds and 
the Milky Way in the parametrization of Pei| ( [1992[ ). 

All synchrotron models yield acceptable fits to the data, and 
there is mild degeneracy between the parameters, in particular 
in the case of a broken power-law continuum with a smooth 
turnover from low to high-energy spectral index. The statisti- 
cally preferred fit (combined statistics^lof 302 for 367 d.o.f.) 
is obtained, however, with a single power-law extending from 
the NIR to the X-ray band. In the case of a broken power- 
law continuum the improvement in the fit is small (combined 
statistics of 299 for 366 d.o.f), and the best-fit break energy is 
located at around 3 keV. With respect to the derived value of 
Ay, the single and broken power-law fits are thus equivalent. 



abundances from |Anders & Grevesse| l 1989) here also for 
3er for a direct comparison with values published in the 



28 We use solar ; 
the intrinsic absorber i 
literature. 

29 We simultaneously use Cash-statistics for the X-ray data and x 2 - 
statistics for the optical/NIR measurements. 



The SED fit with a single power-law is also consistent with 
the similar temporal evolution of the optical/NIR and X-ray 
afterglow light curves (see Figure Al I, This provides strong 



support to the conclusion that both optical/NIR and X-ray data 
probe a common spectral range of the afterglow's synchrotron 
radiation. 

The resulting NIR-to-X-ray SED with the single power-law 
model is shown in Figure [T] Given the modest amount of 
reddening and inferred visual extinction Ay, different extinc- 
tion laws give comparable results. The best-fit is obtained 
in an SMC parametrization (y 2 = 2.9 for 5 d.o.f. in the 
relevant wavelength range), and there is no evidence for the 
presence of a 2175 A dust feature (x 2 = 7.2 for LMC, and 
X 1 - 16.4 for MW extinction curves). The lack/weakness 
of a 2175 A dust feature is corroborated by the X-shooter 
spectrum. While the spectral normalization and large-scale 
shape is set by our (model-dependent) flux calibration proce- 
dure, the inferred slit-loss is uniformly continuous over the 
full spectral range and does not obliterate spectral features of 
the size of the 2175 A bump (FWHM ~ 2000 A observed 
frame). When fitting the spectrum usi ng an extinction law 
in the analytical form based on Fitzpatrick & Massa (1990, 
2007), there is again no evidence for a 2175 A dust feature. 
The C3 parameter of the Fitzpatrick & Massa ( 1990) extinc- 
tion curve (representing the strength of the 2175 A bump) is 
zero within errors. This result is similar to that derived from 
a sample of GRB afterglows with modest ((Ay) ~ 0.4 mag) 
dust obscuration (Schady et al. 2012) and consistent with the 
absence of diffuse interstellar bands (DiBs^j) in our spectrum 
(e.g., |Cox et al.|2006||Weity et al.|2006] ). 

Our data do not constrain the total-to-selective reddening 
Ry, and we report values based on the SMC value (Ry - 
2.93). Best-fit parameters in the single power-law, SMC case 
are A v = 0.15 ± 0.02 mag, N(H) X = (3.5 ± 0.2) X 10 21 cnT 2 
at z - 2.36 and a spectral index of B = 0.78 ± 0.01, which are 
very typical values for GRB afterglow s in general (e.g., Kann 
|et al.|2006|[20T0l|Greiner et al.|201 1) . Also the ratio between 
X-ray absorption and SED-derived extinction is fully within 
the range typically measured for Swif t GRBs (e.g., [S chady 
[eTaLpOO?, 2010; Nardini et alJ|20T0] jWatsonet al. 2013). 
LMC and MW models yield comparable values within errors 
for all parameters. 

We caution and stress that these values are strictly upper 
limits to the host's visual extinction and X-ray absorbing col- 
umn, becaus e of th e presence of strong foreground absorbers 
(see Section [3.2. 1| ), which contribute to both Ay and A^(H)x. 
The expected reddening introduced by the very strong Mg II 
absorber at z = 1.539, for example, could alone be sufficient 



to explain the observed reddening in the data (e.g., Menard 
et al.|2008). We will thus conservatively consider the derived 



values for Ay and N(H)x from the broad-band SED as upper 
limits on the respective parameters in the GRB system in the 
following. 

3.2. Optical/Near-Infrared Spectrum 
The X-shooter spectrum of the afterglow of GRB 1208 15 A 



(Figures Bl B2 and B3( displays a multitude of absorption 



lines. They are interpreted as Lya, LyB, different neutral, low 
and high-ionization metal lines, several fine-structure transi- 
tions, as well as absorption by molecular hydrogen includ- 

30 The /16284 DIB is the only transition of the most prominent DIBs at 
7, = 2.36 that is not located in a window of strong telluric absorption. We set 
a limit on its equivalent width of W t (A 6284) < 0.6 A. 




•< 



Observed wavelength (A) 

Figure 2. X-shooter spectrum around the damped Lya absorption. Grey is 
the X-shooter data and noise level, both binned for clarity. The red-shaded 
region is the f inal de scription and uncertainty of the DLA, using the HI-model 



from Section 



3.2.7 



and a column density of log(iV(Hi)/cm~ ) = 21.95 ± 
0.10. The black dashed lines show a DLA with the same column density, but 
without the vibrationally-excited transitions of molecular hydrogen. 

ing its vibrationally-excited transitions. The GRB redshift de- 
rived from these absorption lines, z — 2.358, was previously 
reported by Malesani et al. (2012). The hydrogen column den- 
sity is log(JV(Hi)/cnr 7 T= 2T55 + 0.10 ( Figure |2J. Here we 
used the model of W 2 (see Section 3.2.7 i, to estimate its in- 
fluence on the DLA strength. In addition, emission from the 
[O m](/l5007) transition is also detected in the NIR arm of 



the spectrum (see Section 3.2.2 and Figure[3]l 



3.2.1. Intervening Systems 

We identify five different intervening systems in the X- 
shooter spectrum that are characterized by at least one tran- 
sition with W b s > 1 A. In detail, there is a very strong Mg II 
absorber at z = 1.539 that also shows lines from various other 
metal transitions (e.g, Fe II, Mg I, Al III). Its rest-frame equiv- 
alent width of W t (A 2796) = 6.1 ±0.2 A makes it the strongest 
intervening Mg II ab sorber yet observe d towards a GRB (e.g., 
Ellison et al. 2006; Verga niet al.|2009 1. In addition, there are 
weaker systems at z = 1.693 (W T (A2 '96) = 0.50 ± 0.03 A) 
and z = 2.013 (W r (A2196) = 0.48 ± 0.04 A). Absorption 
from the C IV doublet is present at z - 2.000 (W r (A 1548) = 
0.63 ± 0.04 A) and z = 2.338 (W r (A 1548) = 0.43 ± 0.04 A). 
The intervening systems will not be discussed further in this 
work. 

3.2.2. Emission Lines 

We used spectral point spread function subtraction (e.g. 
|M0ller|l2OOOl |Fynbo et al.l[20T0l > of the afterglow light to 
search for rest-frame optical emission lines in the X-shooter 
data. In the afterglow-subtracted spectrum (see Figure 13), 
we detect [O III] {A 5007) in emission with a significance of 
6cr. None of the emission lines of Lya, Ha, H/3, [O II], 
[O m](/l4959), [N II] or He 110*4686) are robustly de- 
tected. The Gaussian centroid of [O III](/15007) provides 
a systemic redshift of z - 2.3586 ± 0.0001 for the galaxy 
hosting GRB 1208 15 A. The line width is relatively well re- 
solved, and from the observed full-width-half-maximum of 
FWHM ~ 70kms _I , we derive a velocity dispersion cr Ye \ = 
21^gkms~ 1 after correcting for the instrumental resolution, 
and taking into account its uncertaintjpH Lacking a mea- 
surement of the visible extent of the GRBnost, no dynamical 

3 'Allowing values between the assu med resolving power and its lower 
limit derived from arc-lamps, see Section|2.3| 
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Figure 3. One- and two-dimensional, afterglow-subtracted X-shooter spec- 
trum centered around the location of [O III](/15007) at z = 2.3586 (i.e., ob- 
served ~ 16800 A). The one-dimensional spectrum (black data) fitted by a 
Gaussian (black line) is plotted against the y-axis, while the scaling of the 
two-dimensional spectrum is shown by the color-bar on the right side of the 
plot. The scale along the spatial axis is 072 per pixel. The two-dimensional 
spectrum has been smoothed by a Gaussian kernel with a FWHM of 3 pixel 
(1.5 A) along the dispersion x-axis, and 2 pixel (i.e., 074) along the spatial 
y-axis for illustration purposes. The grey-shaded area represents the spectral 
resolution. 

mass can be derived here. A comparison to GRB hosts (e.g., 
Savaglio et al. 2009, Kriihler et al. 2012a) or len sed galaxies 



at comparable redshift (e.g., Richard et al. 201 1 , Christensen 
et al. 2012), however, indicates a dynamical or luminous mass 



well below 10 lu M o for the host of GRB 1208 15A. The to- 
tal flux (after matching the NIR spectrum to the afterglow 
model, and including its uncertainty) in the line is F([0 III]) = 
(1.8 ± 0.4) x 10~ 17 ergcirT 2 s -1 . The 3cr-confidence limit on 
the Ha flux is F(Ha) < 4.2 x 10~ 17 erg cirT 2 s -1 , correspond- 
ing to a limit on the observed, i.e., not extinction corrected, 
star-format ion rate S FR < 9 M Q yr" 1 u s ing th e conversion 
from |Kennicutt| ( |1998| l with a |Chabrier| ( |2003| > initial mass 
function. 

3.2.3. Absorption Lines, Velocity Structure and Abundances 

The absorption lines in the X-shooter afterglow spectrum 
are relatively narrow, spanning AV = 83^ 2 ° k m s" 1 based on 
Si II(/tl808) and following the definition in Ledoux et al. 



(2006). The absorption is concentrated into two components 
separated by ~ 27 km s _1 and, within our resolution and S/N, 
no further absorption component is required to adequately 
fit the individual absorption lines (Figure R}. The redshifts 
of the blue and red component are z - 2.3579 ± 0.0001 
and 2.3582 ± 0.0001 respectively, which is 63 kms^and 
36 km s _I blueshifted with respect to the systemic redshif t de- 
fined by the [O III](/15007) emission line (Section [3?2?2 i. 

To derive column densities, we fitted Voigt profiles con- 
volved with the instrumental resolution to unsaturated ab- 
sorption lines within VPfit V.10.6P 2 ] Here, we linked 
the broadening parameter b between the individual fits of 
weak absorption lines, and allowed for a small offset (< 
5 - lOkms 1 ) in the wavelength calibration between line fits 
of individual transitions of the same ion. We further assumed 
that only turbulent broadening contributes to b. Finally, we 
imposed the same velocity structure of two central compo- 
nents for all lines, independent of their excitation level or ion- 
ization state. While this is supposedly a simplification of the 
ISM velocity structure, our data do not yield more detailed 
constraints (Figure HI. 

The broadening parameters are constrained to b\ - 
9.7 ± l.Okms -1 for the primary, red component, and bt - 

32 http : //www . ast . cam . ac . uk/~r f c/vpf it . html 
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Table 1 

Column Densities and Abundances 
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Figure 4. Absorption lines and Voigt-profile fits in the GRB 120815A af- 
terglow spectrum. Dark-grey circles and lines show the normalized spectrum 
with statistical errors and the line fit, respectively. Light-grey circles show the 
residuals between data and best fit. The instrumental resolution of X-shooter 
and its error is indicated by shaded, blue areas in panels 2 (VIS arm) and 7 
(UVB arm). The position of individual metal lines of the two absorbing sys- 
tems is indicated by vertical dotted lines. Absorption by strong telluric lines 
is marked with a circled plus (panel 3). 
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Fell* 


2333,2612 
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NiII 4 F 9/2 


2217,2316 


13.23 ± 0.05 




CrII 


2056, 2062, 2066 


13.75 ± 0.05 
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2026 


13.54 ± 0.05 




CI 


1560, 1656 


13.41 ±0.11 
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1306 


15.29 ±0.18 
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1238, 1242 


:£ 14.80 ±0.21 





Note. — ™ For H 2 and CO, this column denotes the different rotational levels of 
the vibration ground-state in the case of H 2 , or the different CO bandheads, respec- 
tively. 

* b ' Total column density as the sum of both velocity components (see Section |3.2.3[ 
and FigureRl. 

* c * Assuming a negligible ionization correction, i.e. A'(Fe) = N(Fe II), N(Zn) 
= N(Zn II), etc., and using the total column density of hydrogen as N(H) = 
N(Ui) + 2N(H 2 ). 

* d ' Summed over all detected states of the same ion. 
^ e ' Quoted uncertainties are the formal errors provided by VPfit. 



11.9 + 3. 5 km s~' for the weaker, blue component, where 
the quoted uncertainties on b are the formal errors provided 
by VPfit. At the given S/N and instrumental resolution of 
FWHM = 28 kms~' in X-shooter's VIS arm, features with 
b > 6kms~ 1 (and thus the velocity profiles of the metal ab- 
sorption) are resolved. The broadening parameters are close 
to the instrumental resolution in the UVB arm (Figure HI). 
The simultaneous fit of all weak lines (excluding the mildly- 
saturated Si ll(A 1808)) has a^ 2 = 1058 for 954 degrees of 
freedom. Si ll(A 1808) is fitted after fixing its line properties 
to the previously-obtained values from weaker lines. 

Given the limited resolution and small separation, the ve- 
locity components are blended. The column densities of the 
single components are correlated and accordingly have sub- 
stantial errors individually. The total column as the sum of 
both components, however, is typically well constrained, and 
we report only total column densities in the following. Our 
best constraint on the metallicity of the GRB-DLA comes 
from the abundance of zinc, an element generally consid- 
ered as being very little affected by dust depletion. Based 
on the well-detected Zn ll(AA 2026, 2062) transitions, we de- 
rive [Zn/H] = -1.15 + 0.11. Estimates from Si II and S II 
agree well within errors. Velocity profiles of selected lines 
are shown in Figure HI and column densities are given in Ta- 
ble m 

All metal transitions bluewards of A KSt ~ 1700 A are 



blended with H* 2 lines (Section 3.2.7 1. Consequently, we con- 
servatively consider the values derived for weak lines that are 
in regions of strong HJ transitions (in particular S ll(A 1250) 



and to a lesser extent also N V(AA 1238, 1242)) as upper lim- 
its only. In addition, absorption lines bluewards of A iest < 
1260 A are located in the damping wing of Lya, adding fur- 
ther to the error of the column density because of the uncer- 
tainty in their continuum placement. 

3.2.4. Depletion of Refractory Elements 

The measured column densities and abundances (Table [TJ 
of the refractory elements such as Fe, Cr, Ni or Mn with re- 
spect to Zn, Si or S indicate that they are depleted onto dust 
grains. For GRB 120815A, we derive a relative abundance 
of [Zn/Fe] = 1.01 + 0.09, a column density of iron locked- 
up in dust of log(JV(Fe) dus 7cirT 2 ) = 16.26 + 0.08 (following 



Vladilo et a l. 2006), and a dust-to-gas ratio based on zinc 
'* Zn = 10l^/ M J(l-10 [Fe/Zn] ) of log(* Zn ) = -1.19 ± 0.12, 
approximately 40% of the value of the SMC. The column 
density of iron in dust grains is a factor of 40 larger than 
the threshold of log(Af(Fe) dust /cnr 2 ) > 14.7 above which 
H2-molecules ar e typically found in QSO-DLAs (e.g., Noter- 
Idaemeetal . 2008). 

Using the SMC as a benchmark^] the depletion-derived ex- 
tinction is Ay - 0.12 ± 0.03 mag. The metals -to-dust ratio 
(log(JV(Hi)/cirT 2 ) + [Zn/H])/A v = (5±l)xl0 21 cm^mag" 1 
is roughly 2.5 times the value of the Local. Different deple- 
tion methods (see, e.g., Savaglio et al. 2003) and comparison 
values from the Local Group allow visual extinctions up to 
Ay ~ 0.6 mag, and consequently a metals-to-dust ratio based 
on the depletio n analysis that is consistent with the G alactic 
value (see e.g, |Bohlin et al.|1978||Predehl & Schmitt|1995| ). 



3.2.5. Molecular Hydrogen 

The most salient feature of the X-shooter spectrum blue- 
ward the Lya absorption is the presence of periodic regions of 
negligible flux (see FigureB] top panel). They are interpreted 
and well described by the Lyman- Werner absorption bands of 
molecular hydrogen, rotationally-excited states between the 
Lyman-limit and a rest-frame wavelength of around 1 1 10 A. 

The fitting and analysis of the molec ular hydrogen transi- 
tions follows Ledoux et al. ( 2002 2003 ) and is shown in Fig- 
ure [3] Given the available data, we can constrain the col- 
umn densities of the rotational levels J — up to J — 3 
in the first five Lyman bands, L - 0-0 to 0-4, which 
yield a total H 2 column of log(JV(H 2 )/cirT 2 ) = 20.53 ± 0.04. 
Details of individual rotational transitions are provided in 
Table [T] Together with the column density of neutral hy- 
drogen log(A^H/cm -2 ) = 21.95 ± 0.10, this translates into 
a molecular fraction of log/ = —1.14 ± 0.10, where f = 
2N(U 2 )/(N(Ui) + 2N(H 2 )). The kinetic temperature for the 
excitation of the J — 1 state is T\^ n = 200^ ° K, which is 



higher than observed locally (e.g., Savage et al. 1977[|Turnliri-| 

fn et al. 2002 ), but in the range of QSO-derived temperatures 
rianand et al'|2U05) . 

The best-fit redshift and broadening parameter of the 
molecular hydrogen lines are z(H 2 ) = 2.3582 and b(H 2 ) = 
8.7 ± 0.6kms~' , consistent with the pr imary component of 
the metal-line absorption (Section |3.2.3| l. 

3.2.6. Lack of Ca rbon Monoxide 



This requires assumptions on the SMC's dust depletion ([Fe/Zn] S MC = 
- 1 .0), and metallicity ([M/H] SMC = -0.7), the average reddening per N(H i) 
column, and the total-t o-selective reddening, which were all taken from 
IProchaska et al.|<2007al and references therein. 



We do not detect absorption signatures from the CO 
molecule. While there is an absorption line at ~ 4960 A with 
observed equivalent width W b s = 0.64 ± 0.10 A, conspic- 
uously close to the CO AX(2-0) bandhead at z = 2.36, we 
consider it unrelated to CO, as none of the other expected 
CO transitions are detected. Also, the CO redshift would be 
significantly offset from the value of either absorption, emis- 
sion or H 2 . We thus conclude that this single line is not 
caused by CO, but related to an intervening absorber. Given 
the presence of a second absorption line at ~ 4966 A with 
W obs (A 1550) = 0.33 + 0.07, the most likely interpretation of 
these features is the C IV doublet at z — 2.2042. 

Using the wavelength range of the six strongest CO AX 
bandheads (CO AX(0 - 0) to CO AX(5 - 0)) redshifted to z = 
2.3582 as reference, we set a 3<x upper limit on the total CO 
column densityp]of log(JV(CO)/cirT 2 ) < 15.0. This value is 
fairly independent on the CO excitation temperature (assumed 
to be greater than Tqmb — 9 K at z - 2.36) and CO broadening 
parameter, because the latter would not be resolved for typical 



values with our X-shooter data (e.g., Noterdaeme et al. 2010| ). 

The N(CO) to N(U 2 ) ratio is smaller than 10 3j , which 
is typical for Galactic measurements with / < 0.1 (Burgh 
et al. 2010). The H 2 -bearing system towards GRB 120815A 
is thus characterized as a diffuse cloud, consistent with the 
small amount of re ddening and extinctio n (in comparison to 
Galactic sightlines; [Snow & McCall|2006|). 

Because of the sub-solar GRB-DLA metallicity, 
N(CO)/N(H 2 ) could even be below what is typically 
observed in the Milky Way. The CO luminosity per unit 
luminosity or per unit SFR of local dwarf galaxies, for exam- 
ple, is significa ntly lower tha n extrapolated from metal-rich 
spirals (Schruba et al. 2012). This can be interpreted as 
evidence for a strong metallicity dependence of the CO to 
H 2 ratio below half of Z (see also, e.g., |Leroy et al.|2011[ ). 
Measurements from GRB-DLAs similar to the one presented 
in this work will help to probe and study the CO abundance 
and the H 2 to CO ratio over a broad range of metallicities (or 
stellar masses) and redsh ifts directly and co mplementary to 
conventional studies (e.g., Genzel et al.|2012 1. 

Our upper limit on N(CO)/N{H 2 ) is comparable to or lower 
t han derived from C O detections in high-redshift s ystems 
(|Srianand et al.|2008[rProchaska et"al . 2009 ; Noterdaeme et ah] 



2010| i, which we readily relate to the lower metallicity and 
dust-extinction of the GRB 1208 15 A DLA as compared to 
the sightlines in which CO was seen. In QSO-DLAs without 
CO detection, the constraints o n CO/H? extend dow n to limits 
of N(CO)/N(H 2 ) < 10" (e.g.,|Petitjeanetal.|2002|). 



3.2.7. Detection of Vibrationally-excited H 2 

After establishing the presence of molecular hydrogen, we 
search for vibrationally-excited H 2 by usin g the synthetic 
spectru rrp] calculated in Draine (2000) and |Draine & Hao 
(2002). We downgraded a theoretical R — 10000 spectrum 
to the resolution of X-shooter's UVB arm (R ~ 6000), limited 
the wavelength range to regions above the red damping wing 
of Lya, and excluded strong metal absorption lines caused 
by the GRB-DLA and the intervening systems in the further 
analysis. 

Cross-matching the downgraded spectrum with our data 
(Figure [6]l using different optical depths with respect to the 
input model and varying redshifts, we find a single and well- 



Morton & Noreau 



34 Using molecular data for CO from 

35 http : //www . astro . pr inceton . edu/~draine/grbh2 .html 
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Figure 5. X-shooter spectrum between 3300 and 4200 A illustrating the presence of lb absorption. The uppermost, blue-framed panel shows an overview, 
while the lower 5 panels show zoom-ins of 100 A each below 3800 A, where most of the Lyman-Wemer absorption bands are located. Light-grey lines are 
the normalized spectrum, while dark-grey lines indicate the error spectrum. Red lines denote the best-fit l b. model, with i ndividual transitions marked with red 
labels. In the top panel, we also mark prominent metal absorption lines previously detected in GRB-DLAs iChristensen et al. 201 1 1 and those of the intervening 
systems. In the lower panels, JO, Jl, and so forth denote transitions from the J = 0, J = 1 rotational level of the vibrational ground state of the Lyman-Werner 
bands of H2 . 
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Figure 6. X-shooter spectrum between 4100 and 5600 A showing the presence of Hj lines. The uppermost, blue-framed panel shows an overview over the whole 
spectral range, while the lower 5 panels zoom in of 1 80 A each. We selected regions where most of the absorption bands of HI, are located. Light-grey lines 
always show the normalized spectrum, while dark-grey lines indicate the error spectrum. Red lines denote the best-fit H*, model. In the top panel, we also mark 
prominent metal absorption lines previously detected in GRB-DLAs ( Christensen et al. 201 1 1. In the lower panels, several individual H^ transitions are identified 
using standard nomenclature with lower and upper vibrational and rotational quantum numbers. 
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the mod el, we use the optical afterglow light curve (see Sec- 
tion 2.2 1, the after glo w spec tral slope jS, the inferred extinction 
(Sections 3.1 and|3.2.4[) and the broadening parameter b of 
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Figure 7. Logarithmic Ax distribution for different H^ model redshifts and 
HJJ absorption line strengths (the ratio between the applied optical depth and 
the one in the input model). 

defined minimum in the^- 2 -space at a redshift consistent with 
the ground-state H2 and metal transitions (Figure[7ji. The best 
match is obtained at z — 2.3581 ± 0.0001, and the presence of 
H*, as compared to an unmodified spectrum is unambiguous 
(FigurelTli. 

While the overall match between the input spectrum and 
data is reasonable, there are subtle differences between the 
data and synthetic model. These are evident in the different 
line ratios between vibrational levels (see, e.g., the expected 
but undetected strong line at ~ 4560 A in Figure 6j). This is 
likely related to different initial conditions — suchas parti- 
cle densities, cloud-afterglow distance (see Section |4~T) , dust 
shielding, afterglow luminosity and spectrum, and/or back- 
ground UV-radiation field — with respect to the assumptions 
in the theoretical calculation or uncertain transition probabili- 
ties between different excited levels of the H2 molecule. In 
particular, we derive a distance between the absorbin g ga s 
cloud and GRB 1208 15A of d ~ 0.5 kpc (see Section pUj, 
much larger than the d ~ 1 pc assumed in the model. This 
will inevitably lead to different line strengths and ratios be- 
tween the various H*, transitions. 

4. DISCUSSION 
4.1. Distance between GRB and Absorbing Cloud 

Our observations indicate that the metal-lines and H2 ab- 
sorption are not caused by circumburst material. The presence 
of H2 alone sets a lower limit on the burst-cloud distance o f at 
least d > 10 pc ( |Whalen et al.|2008HLedoux et al.|2009l >. In 
addition, the detection of neutral species at the same redshift, 
in particular Mg I (Figure 0), is generally interpreted as an 
indication that t he gas is at substantial distances (d > 50 pc) 
from the GRB (Prochaska et al.|2006| ). 

These initial considerations are confirmed by modeling 
the photo-excitation of the Fe II*, Fe II**, Fe II 4 F 9 /2 and 
Ni n 4 F9/2 transitions^ Given the small separation in veloc- 
ity of both absorption components, and the limited resolution 
of our X-shooter data, we onl y model the tota l column densi- 
ties (see also the discussion in Vreeswijk et al.|20"T3 1. 

The photo-excitation modeling follows closely the method- 
olo gy presented an d applied in |Vreeswijk et al. ( |2007| [20T3] l 
and |Ledoux et al.| ([2009), and we refer to these publications 
for a detailed description of the technique. As an input to 

36 In the modeling, we excluded Si II*, because excited states of Si II are 
also seen in very actively star-forming galaxies (e.g., |Pettini et al.|2002) . Si II 
is easier to excite than Fe II or Ni II, and a fraction of the Si II excitation 
might thus not be due to the GRB afterglow radiation. The observed ratio of 
N(Si II*)/W(Si II) ~ 1% could, for example, be provided by electron densities 
around lOcirT 3 in a T ~ 10 4 K medium, or UV-backgrounds about 10 3 times 
that of the Galaxy fSilva & Viegas|2002} . 



the absorption lines ( |3.2.3[ ). The dust is placed inside the ab- 
sorbing cloud that is being excited, and the distance d from 
the GRB is defined to the near-side of the cloud. 

A good fit (y 2 /d.o.f = 2.3/4) is obtained for a distance be- 
tween gas cloud and burst of d — 4 60 ± 60 pc in the case of 
atomic data from CLOUDY ( |Verner e t al. 1999). A similar re- 
sult is obtained (y 2 /d.o.f = 4.3/4, d - 530±80pc) when using 
atomic data as described in Vreeswijk et al. (2007). For a de- 
tailed comparison between the different sets of ato mic data in 
this context, we refer to the discussion in Ledoux e t al.| ( [2009l >. 
Given the lack of a spectral time series, and relatively few 
measurements of excited metal transitions, the cloud size is 
unconstrained by our data. The neutral gas seen in absorption 
in the GRB-DLA is thus not directly related to the GRB pro- 
genitor's immediate environment or stellar wind-driven mass 
loss. 

4.2. Detection of High-Ionization Species 

We also detect the imprint of the highly-ionized transition 
N V (ionization potential for creation of 77.5 eV) with a high 
column density at a redshift consistent with the other metal- 
lines. N_Vis a common feature observed towards GRBs ( |Fox| 
et al. 2008) and is suggestive of a circumburst origin in the 
vicinity of the GRB (Prochas ka et al.||2008a[ ). N V is, how- 
ever, also present in UV spectra of O-type stars, galaxies and 
QSO-DLAs (e.g., IPetfini et al.|2002||Fox et al.|2009) . In the 
case of GRB 1208 15 A, in particular, the N V measurement 
is, due to blending with H* 2 , strictly an upper limit, suffers 
from low S/N and has thus large errors on its column density 
as well as uncertainties in its velocity structure. Our detec- 
tion of N V is broadly consistent with the measurements for 
other GRB-DLAs, and does not necessarily point to an origin 
in the circumburst environment, but could also be of interstel- 
lar origin as argued in |Fox et al. (2008). A location in the 
ISM is supported by the non- variability of N V in a spectral 
time-series for GRB 0803 10 ( |De Cia et r al.pOT2] >. 

4.3. Molecular Hydrogen in GRB-DLAs 
4.3. 1 . Metallicity and Neutral Hydrogen Column Density 

The column densities of H 1, H2 and the molecular fraction 
/ of the DLA towards GRB 1208 15 A are high when com- 
pared to the full sample of QSO-DLAsF] The high JV(Hi) 
is, however, very typical for GRB-DLAs. When consider- 
ing the significant dust depletion and metal abundance along 
the GRB 1208 15 A sightline, a consistent picture is emerging: 
High molecular fractions (log(/) > -3) are exclusively seen 
at significant metallicities ([M/H] > —1.5), while they are not 
restricted to the very high-end of hyd rogen column densi ties 



Ledoux et al. 



2003) 



of log(JV(Hi)/cirr 2 ) ~ 21.0 (see, e.g.. 

The high column densities of neutral hydrogen typical for 
GRB-DLAs are thus not the primary condition for the pres- 
ence of H2 in DLAs. This is illustrated in Figure [8] which 
shows the incidence of H2-bearing gas as functions of Hi 
column densities and metallicity. While there is an obvi- 
ous discrepancy in the distributions of JV(Hi) between GRBs 
and QSO-DLAs, the difference between H2-poor and H2-rich 

37 Noteworthy exceptions of dusty and metal-enriched QSO-DLAs (e.g., 
Fynbo et al.|201l) or QSO-DLAs with proper ties similar to those of typical 



GRB-DLAs exist (e.g., Guimaraes et 



ith properti 
t al.|2012) . 
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Figure 8. The central panel shows metallicity versus hydrogen column den- 
sity for QSO and GRB DLAs in which H2 was extensively searched for. 
Circles are QSO-DLAs, where red and yellow circles correspond to DLAs 
with and without significant molecular hydrogen. Blue data are GRB-DLAs 
and the blue stars are GRB 120815A and GRB 080607, where H 2 has been 
detected. The grey background shading indicates a region where many GRB- 
DLAs are expected, but are not studied at high spectral resolution. The upper 
and right panels illustrate the respective normalized histograms. Blue hatched 
histogram are GRBs including both H2-poor and ILj-rich sightlines due to 
low number statistics. Red-filled and yellow histograms are QSO-DLAs with 
and without molecular hydrogen, respectively). Data for QSOs are from No- 
|terdaeme et al.|(2008 i; Fynbo et al. (2010 2011 1; Guimarae Fet al.|<2012>, 
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Figure 9. Depletion of QSO and GRB DLAs versus their metallicity. Sym- 
bols, colors and shadings are the same as in Figure [8] GRB 080607 is not 
shown in this plot beca use no constraints on the depletion factor are available 
(Prochaska et al.|2009) . 

QSO sightlines with respect to their N(H 1) is much more sub- 
tle (Figure [8] top panel). 

4.3.2. The Role of Dust 

Dust plays a major role for the presence of molecular gas by 
catalyzing the formation of molecules on the surface of dust 
grains as well as shielding against Lyman-Werner photons. 
For GRB 1208 15A, the presence of dust is evident from the 
relative abundance of Fe, Ni, Mn and Cr with respect to Zn 
([Zn/Fe] = 1.01 + 0.09). This value is larger than the relative 
abundance of previous GRB-DLAs observed at high spectral 
resolu tion, for w hich elaborate H2-searches were performed 
( |Ledoux et al.|2009) >. 



Similar to the case of metallicity, there is a higher proba- 
bility of finding molecular hydrogen for more dust-depleted 
QSO sightlines (Figure [9]). There seems to be a lower limit 
of [M/Fe] of ~ 0.3 below which no strong signatures of 
H2-molecules (neither towards QSO nor GRBs) are observed 
(Figure 15). The values of metallicity, dust depletion and 
N(H 1) are effectively combined into a density of iron locked 
up in dust Af(Fe) dust , which has been ide ntified as a pri- 
mary driver of H2 detection in QSO-DLAs (Noterdaeme et al. 
2008). 



Guimaraes et al. (2012 1 study the QSO-selected system 
J0816 and identify it as an analog to GRB-DLAs with very 
similar N(Hi) and [Zn/H] as measured for the DLA of 
GRB 1208 15A (Figure |8). The column density of N(H 2 ) 
in the QSO-DLA, and thus the molecular fraction /, how- 
ever, is two orders of magnitude lower than derived for 
GRB 1208 15A. This could possibly be related to the lower 
depletion. The mean relative abundanca^j in the QSO-DLA 
was [Zn/Fe] = 0.48 + 0.06, while we observe [Zn/Fe] = 
1.01 ±0.09 (Figure [5J. 

4.3.3. Implications for the Population of GRB-DLAs 

The region of elevated metallicity, depletion and N(H 1) is 
exactly the parameter space where we expect many of the 
more dust-extinguished, and thus UV-fainter afterglows to be 
located. This combination directly leads to a high absolute 
Af(Fe) dust , and thus an efficien t conversion from H 1 to H2 on 
the dust grain surface (Schaye 2001 Cazaux & Spaans 2004). 
While a high column density of dust is a necessary condl- 
tion for creating molecular hydrogen on dust grains, the ob- 
served H2 column density is ultimately set by the interplay 
between H2 formation and photo-dissociation, for example 
by the background UV-radiation field in the GRB host. A 
lack of molecular gas in GRB-DLAs with respect to QSO- 
DLAs of similar ^V(Fe) dust would thus indicate that photo- 
dissociation of H2 is more prominent in GRB hosts than it is in 
the galaxy counterparts of GRB-DLAs. In such a scenario, Fb 
needs to be absent already prior to the burst, because the de- 
rived distances between the absorbing gas and GRB are much 
larger than the region of d < lOpc, in which the afterglow's 
UV-flux photo-dissociates H2. The previous lack of molecu- 
lar gas in GRB-DLAs, however, is co nsistent with their ob - 
served metallicity and dust depletion ( Ledoux et al.|[2005) , 
and does not necessarily point to a difference between GRB 
and QSO-DLAs with respect to the presence of H2. In fact, 
GRB 050820A is still the only H 2 -less GRB-DLA with phys- 
ical properties that are similar to those of H2-bearing QSO- 
selected systems (Figure 15). 

The measured hydrogen column density, metallicity and 
dust-depletion of GRB 1208 15 A are fairly representative of 
GRB -DLAs prev iously observed at lower spectral resolution 
(e.g.,|Savaglio et al.|2003||Savaglio 2006]|Fynbo et al.|2009 



|Rau et al.||2010| |Christensen et al.||20li) . Based on this 
low-resolution data, we expect the distributions of the re- 
spective physical quantity to be centered roughl)p1 around 
log(JV(Hi)/cnr 2 ) ~ 21.6, [M/H] ~ -1 and [M/Fe] ~ 0.8. 
This parameter space is illustrated by the grey shaded area in 
Figures [8] and I5J Many GRB-DLAs have physical parameters 
that are favorable for H2 to exist, but their molecular con- 
tent could not be probed with afterglow spectra of low S/N, 



Converted into our metallicity scale based on |Asplund et al. J2009 '. 
39 We note that also in the low-resolution sample, biases towards lower 
metallicity and lower depletion values are present I Fynbo et al. 2009 1. 
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low resolution, or both. Once we have gained better observa- 
tional access to similar afterglows through detailed observa- 
tions with sensitive spectrographs of high-enough resolution 
(such as X-shooter or UVES), it seems likely that a higher 
fraction of those GRBs will also show the presence of molec- 
ular gas. 

4.4. Vib rationally-excited H2 in GRB-DLAs 

The presence of vibrationally-excited states of molecular 
hydrogen (H*) along GRB sig htlines was first postulated by 
|Draine| ( |2000| ) and detailed in |Draine & Hao] ((2002). The 
strong features, however, that were expected to be produced 
by the vibrationally-excited H2 from the GRB's birth cloud, 
and detectable in even low-resolution spectra, were never 
clearly observed. Excited states would nevertheless be present 
in cases where the m olecular cloud only intersects the GRB 
sightline (Draine & Hao 2002 ) at small enough d istances. In- 
deed, while in the work of [Draine & Hao| ( |2002p the total H* 
column density produced inrough UV-pumping is N(HX) - 
10 19 cirT 2 , the best overall match is obtained here with a fac- 
tor 6 decreased optical depth of the Hi transition s compared 
to the input model. In the case of GRB 080607, [Shefferet al. 
(2009) perform a dynamic modeling of the H* 7 absorption, and 
derive a distance between GRB 080607 and its H2-bearing 
cloud between 230 and 940 pc. This is comparable to our es- 
timate of the distance between GRB 1208 15 A and its DLA of 
0.5±0.1kpc. 

As shown in Figures |2| and [6] the H^ transitions af- 
fect the measurement oflv(Hi). If H2 and W 2 are more 
common in GRB-DLAs as speculated above, JV(Hi) could 
have been overestimated in previous cases. In the case of 
GRB 1208 15 A, the H 2 transitions decrease the best-fit col- 
umn density by ~0.15 dex. Similarly, blending with H* 2 could 
be a serious concern for intrinsica lly weak metal l ines, such 
as S ll(AA 1250, 1253) (see also |Sheffer et al.|2009| l. 

H* 2 opens a route to establish a positive presence of molecu- 
lar hydrogen also in those cases where the observations do not 
cover the wavelength range bluewards of Lya, and thus par- 
ticularly in afterglow spectra of GRBs in the redshift range 
1 < z < 2 . Strong vibrationally-excited levels of H2 extend 
up to /} rest ~ 1600 A, and in contrast to the Ly man-Werner 
bands at A KSt < 1120 A, their identification is not compro- 
mised by the Lyor forest, and could be performed also with 
lower-resolution data. Primary candidates for further H 2 (and 
thus H2) searches are the high-quality, high-resolution spectra 
of low-redshift afterglows (in which the individual lines could 
be identified even if weak), or the spectra of metal- and dust- 
rich GRB absorbers (in which a significant column density of 
H2 could to be present). 

5. CONCLUSION 

We have presented optical-to-NIR X-shooter spectroscopy 
of the afterglow of GRB 1208 15A at z - 2.36, supplemented 
by optical/NIR photometry from GROND and X-ray data 
from the Swift satellite. The bright afterglow emission, ob- 
served through efficient instruments at large telescopes, pro- 
vides a detailed probe of the physical properties of the inter- 
stellar medium in a high-redshift star-forming galaxy. 

The sightline towards GRB 1208 15A is characterized by a 
strong DLA with log(JV(Hi)/cirT 2 ) = 21.95 + 0.10 and sub- 
stantial amount of molecular hydrogen with a molecular frac- 
tion / of 7% (log/(H 2 ) = -1.14 ± 0.10), characteristic of 
Galactic diffuse clouds. This presents only the second unam- 
biguous detection of H 2 in a GRB-DLA, and the first for a 



GRB-DLA with properties very similar to an average GRB 
sightline. In addition, we detect vibrationally-excited states 
of H2, which opens a second route for positive searches of 
molecular gas in GRB afterglow spectra. 

Our measurements of DLA metallicity ([Zn/H] = -1.15 ± 
0.11), relative abundance ([Zn/Fe] = 1.01 ± 0.09) and visual 
extinction (Ay < 0.15 mag) are common among GRB-DLAs, 
and likely typical for the population of GRB-DLAs in general. 
GRB 1208 15 A thus stands in marked contrast to the metal- 
rich, H 2 -bearing DLA of GRB 080607. This illustrates that 
H2 is present in at least a fraction of the average GRB-DLAs 
systems. The detection rate of molecular gas in GRB-DLAs 
could increase once similar afterglows are observed routinely 
at higher spectral resolution. 

Similar to many other physical properties probed by GRBs 
or their hosts, such as the distributions of galaxy brightness, 
mass, star- formation rate (Kriihler'et al.|[201 1 ; Hjorth~et al.| 



2012| Perley et al.||2013), galaxy col or (Rossi et al. _ 

Lya emission ( Milvang- Jensen et al. 2012|l, dust-reddening 



(Cenko et al. 2009; Greiner et al. 2 01 1), soft X-ray absorption 
(e.g., Fynb o et al .120 09 ; Campana et al ]|20T0} JKruhler et al. 
2012b) or metallicity (Prochaska et al. 2009), selection effects 
play an important role when studying H2 in GRB-DLAs. 

A larger sample of afterglows observed in a similar way as 
presented in this work (Fynbo et al. 2013, in preparation) will 
enable further progress in this field and deeper insights with 
respect to the statistical presence and properties of molecules 
along GRB sightlines. Additional detections of molecular gas 
in GRB-DLAs would allow detailed individual as well as sta- 
tistical studies, and, coupled with host follow-up and sub-mm 
spectroscopy, provide unprecedented insights into the process 
and conditions of star-formation at z ~ 2 and above. 
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Table Al 

Optical magnitudes of the afterglow of GRB 1208 15A as measured by GROND 



Time 


Exposure'"' 




Brightness*' 








g'-band 


r'-band 


('-band 


z'-band 


(ks after 7o) (c) 


(s) 


(magAB) 


(mag AB ) 


(magAB) 


(magAB) 


0.169 


35 


18.92 ±0.03 


18.03 ± 0.01 


17.60 ± 0.02 


17.32 ±0.02 


0.273 


35 


18.85 ± 0.03 


17.95 ± 0.01 


17.52 ± 0.02 


17.25 ± 0.02 


0.375 


35 


18.75 ± 0.03 


17.87 ± 0.01 


17.45 ± 0.02 


17.16 ±0.02 


0.482 


35 


18.75 ± 0.02 


17.88 ± 0.02 


17.47 ± 0.02 


17.19 ±0.02 


0.698 


115 


18.85 ± 0.02 


18.01 ±0.01 


17.58 ±0.01 


17.32 ± 0.02 


0.884 


115 


19.02 ± 0.02 


18.16 ±0.01 


17.75 ± 0.02 


17.48 ± 0.02 


1.069 


115 


19.14 ±0.01 


18.27 ± 0.01 


17.85 ± 0.01 


17.60 ± 0.02 


1.256 


115 


19.21 ±0.01 


18.36 ±0.01 


17.95 ± 0.02 


17.66 ± 0.02 


1.451 


115 


19.31 ±0.01 


18.44 ±0.01 


18.04 ±0.02 


17.76 ± 0.02 


1.638 


115 


19.38 ± 0.01 


18.50 ±0.01 


18.08 ± 0.01 


17.80 ± 0.02 


1.834 


115 


19.41 ± 0.01 


18.57 ±0.01 


18.14 ±0.02 


17.86 ± 0.02 


2.027 


115 


19.49 ± 0.01 


18.62 ±0.01 


18.19 ±0.02 


17.92 ± 0.02 


2.361 


375 


19.58 ± 0.01 


18.71 ± 0.01 


18.29 ±0.01 


18.02 ±0.02 


2.807 


375 


19.70 ±0.01 


18.84 ±0.01 


18.38 ± 0.02 


18.13 ±0.02 


3.258 


375 


19.79 ± 0.02 


18.91 ±0.01 


18.47 ± 0.02 


18.18 ±0.02 


3.707 


375 


19.89 ±0.01 


19.00 ± 0.01 


18.59 ±0.03 


18.30 ±0.02 


4.165 


375 


19.95 ± 0.01 


19.09 ± 0.01 


18.65 ± 0.02 


18.38 ± 0.02 


4.607 


375 


20.02 ± 0.02 


19.16 ±0.01 


18.72 ±0.02 


18.44 ±0.02 


5.052 


375 


20.13 ±0.02 


19.21 ± 0.02 


18.81 ±0.02 


18.50 ±0.03 


5.498 


375 


20.18 ±0.03 


19.29 ±0.01 


18.91 ± 0.02 


18.56 ±0.02 


6.037 


4x35 


20.29 ± 0.02 


19.37 ± 0.03 


18.95 ±0.03 


18.67 ± 0.03 


6.465 


4x35 


20.32 ± 0.04 


19.43 ± 0.03 


18.93 ± 0.04 


18.77 ± 0.04 


6.895 


4x35 


20.36 ± 0.05 


19.48 ± 0.02 


19.02 ± 0.03 


18.80 ±0.03 


7.326 


4x35 


20.39 ± 0.03 


19.57 ± 0.02 


19.10 ±0.03 


18.84 ±0.04 


7.763 


4x35 


20.42 ± 0.04 


19.60 ± 0.03 


19.17 ±0.03 


18.97 ± 0.04 


8.196 


4x35 


20.44 ± 0.05 


19.62 ± 0.02 


19.18 ±0.03 


18.93 ± 0.04 


8.631 


4x35 


20.54 ± 0.04 


19.67 ± 0.02 


19.27 ± 0.03 


19.03 ± 0.03 


9.062 


4x35 


20.58 ± 0.04 


19.69 ± 0.02 


19.25 ± 0.03 


19.01 ± 0.03 


9.497 


4x35 


20.63 ± 0.04 


19.75 ± 0.02 


19.27 ± 0.03 


19.13 ± 0.04 


9.927 


4x35 


20.63 ± 0.04 


19.82 ±0.02 


19.37 ±0.03 


19.13 ±0.03 


10.363 


4x35 


20.75 ± 0.03 


19.83 ± 0.03 


19.39 ± 0.03 


19.18 ± 0.04 


10.623 


4x35 


20.79 ± 0.04 


19.86 ±0.02 


19.43 ± 0.03 


19.24 ± 0.04 


11.053 


4x35 


20.81 ± 0.06 


19.91 ± 0.03 


19.48 ± 0.04 


19.22 ± 0.04 



Note. — '"' Integration time of the individual image. Stacked images are given as number of 
images times exposure of an individual image 

' b) All magnitudes in this table are in the AB system and uncorrected for Galactic foreground 
extinction. The quoted error is statistical only, and there is an additional systematic error in the 
absolute photometric calibration, which is estimated to be around 0.04 mag in g'r'i'z' ■ 
(c) To is set as the time of the Swift/BAT trigger, i.e., 2012-08-15 02:13:58 UT. 
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APPENDIX 
A. GROND PHOTOMETRY AND MULTI-COLOR LIGHT CURVES 

The temporal evolution of the optical/NIR photomet ric and X-ray data was fitted simultaneously with a two-fold broken power- 
law, connected smoothly at the break times following Kriihler et al. ( 2009 ). The functional form of this empirical fit is provided for 



example in Schulze et al. ( |2011[ ). The light curve shape is characterized by a shallow, early rise with an index a\ - -0.18 + 0.02, 



where the sign follows the convention that F v (t) <x r a . The optical/NIR afterglow peaks at t\ - 440 + 30 s after the BAT 
trigger time. Afterwards the light curve decays with an index of ai = 0.52 ± 0.01, before another breaking to a steeper decay of 
a-$ - 0.86 ± 0.03 at ti - 4.3*2? ks. The light curve is reasonably well fitted with this phenomenological model (y 2 = 261 for 242 
degrees of freedom), and no chromatic evolution is apparent within our data set. The overall light curve behavior is reminiscent 



of GRB 080710 ( jRriihler et al. 2009}, in particular the achromat ic and shallow, early rise of the light curve in combination with 
;te cted byB./ 
All magnitudes in Figure Al and Tables |A1||A2| are in the AB system and uncorrected for the expected Galactic foreground 



the soft y-ray emission detected 



by BAT (Markwardt et al. 2012). 
and Tables IA11 IA2I are in the AB 



extinction, corresponding to areddening of E B _ V - O.lOmag (Schlafly & Finkbeiner|2011 1. 



B. FULL X-SHOOTER SPECTRA OF THE AFTERGLOW OF GRB 1208 15A 



Table A2 

Near-infrared magnitudes of the afterglow of GRB 1208 15 A as measured by 

GROND 
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Time 


Exposure'"' 




Brightness' 11 ' 










/-band 


//-band 


£ s -band 


(ks after r ) (c) 


(s 




(magAB) 


(magAB) 


(magAB) 


0.191 


6x 


10 


16.84 ± 0.04 


16.59 ± 0.05 


16.28 ± 0.05 


0.295 


6x 


10 


16.80 ± 0.04 


16.54 ± 0.05 


16.28 ± 0.05 


0.397 


6x 


10 


16.73 ± 0.04 


16.48 ± 0.05 


16.19 ±0.05 


0.505 


6x 


10 


16.78 ± 0.04 


16.52 ± 0.05 


16.23 ± 0.06 


0.724 


12 x 


10 


16.95 ± 0.04 


16.75 ± 0.05 


16.50 ± 0.05 


0.909 


12 x 


10 


17.10 ±0.04 


16.85 ± 0.05 


16.66 ± 0.05 


1.094 


12 x 


10 


17.20 ± 0.04 


16.88 ± 0.05 


16.63 ± 0.06 


1.281 


12 x 


10 


17.27 ± 0.04 


16.99 ± 0.05 


16.87 ± 0.06 


1.477 


12 x 


10 


17.36 ± 0.04 


17.21 ± 0.05 


16.86 ± 0.06 


1.663 


12 x 


10 


17.46 ± 0.05 


17.19 ±0.05 


17.01 ± 0.06 


1.859 


12 x 


10 


17.41 ± 0.04 


17.14 ± 0.05 


16.93 ± 0.06 


2.053 


12 x 


10 


17.57 ± 0.04 


17.25 ± 0.05 


17.09 ± 0.06 


2.384 


30 x 


10 


17.60 ± 0.04 


17.45 ± 0.05 


17.09 ± 0.05 


2.830 


30 x 


10 


17.80 ± 0.05 


17.53 ± 0.05 


17.27 ± 0.06 


3.282 


30 x 


10 


17.80 ± 0.05 


17.56 ± 0.05 


17.29 ± 0.07 


3.732 


30 x 


10 


17.90 ± 0.05 


17.69 ± 0.06 


17.57 ± 0.09 


4.188 


30 x 


10 


17.93 ± 0.05 


17.77 ± 0.08 


17.54 ± 0.07 


4.630 


30 x 


10 


18.02 ±0.05 


17.74 ± 0.07 


17.64 ±0.10 


5.076 


30 x 


10 


18.21 ± 0.08 


17.80 ±0.11 


17.63 ± 0.17 


5.521 


30 x 


10 


18.12 ±0.05 


17.98 ± 0.07 


17.82 ±0.10 


6.060 


24 x 


10 


18.21 ±0.07 






6.487 


24 x 


10 


18.28 ±0.10 






6.488 


72 x 


10 




18.07 ± 0.07 


17.90 ± 0.09 


6.917 


24 x 


10 


18.42 ±0.07 






7.349 


24 x 


10 


18.38 ± 0.07 






7.786 


24 x 


10 


18.41 ±0.07 






7.786 


72 x 


10 




18.23 ±0.05 


17.95 ± 0.07 


8.218 


24 x 


10 


18.42 ±0.07 






8.653 


24 x 


10 


18.56 ±0.07 






9.085 


24 x 


10 


18.53 ±0.06 






9.087 


72 x 


10 




18.31 ±0.05 


18.15 ± 0.06 


9.519 


24 x 


10 


18.53 ±0.06 






9.950 


24 x 


10 


18.66 ±0.07 






10.386 


24 x 


10 


18.73 ±0.07 






10.601 


96 x 


10 




18.42 ±0.05 


18.11 ±0.06 


10.823 


24 x 


10 


18.65 ±0.07 






11.251 


24 x 


10 


18.75 ±0.08 







Note. — '"' Integration time of the individual image. Stacked images are 
given as number of images times exposure of an individual image 
* b ' All magnitudes in this table are in the AB system and uncorrected for Galactic 
foreground extinction. The quoted error is statistical only, and there is an addi- 
tional systematic error in the absolute photometric calibration, which is estimated 
to be around 0.06 mag in JH and 0.08 mag in K s . 
(c) T is set as the time of the Swift/BAT trigger, i.e., 2012-08-15 02:13:58 UT. 
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Figure Al. GROND optical/NIR and XRT 0.3-10 keV light curves of the afterglow of GRB 120815A in the top panel. The data were fitted with a three-segment, 
smoothly broken power-law, with residuals shown in the bottom panel. The grey shaded area represents the time interval in which the X-shooter spectroscopy 
was obtained. X-ray data were converted into a flux density at 3 keV using a spectral index of fi = 0.78, and scaled by a factor of 20 to enhance clarity. 
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Figure Bl. X-shooter UVB-arm spectrum of the afterglow of GRB 120815A. Black lines show the spectral data, grey lines the noise level, the green line is the 
GRB afterglo w model, and the red-line shows the DLA modeling. The position of absorption lines that are typically associated with GRB-DLAs (taken from 
Christensen et al.|201 1) are indicated by brown lines and ions. Individual lines associated with intervening absorbers are labeled with colored lines and labels. 
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Figure B2. X-shooter VIS-arm spectrum of the afterglow of GRB 1208 15 A. Lines and labels are the same as in Figure |B 1 1 Additionally, strong telluric lines are 
indicated by circles in the upper part of each panel. 
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Figure B3. X-shooter NIR-arm spectrum of the afterglow of GRB 1208 15 A. Lines and labels are the same as in Figure[B2] 



